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We investigate the possibility ol producing ultracold NH3 molecules by sympathetic cooling in a 
bath of ultracold atoms. We consider the interactions of NH3 with alkali-metal and alkaline-earth 
atoms, and with Xe, using ab initio coupled-cluster calculations. For Rb-NH3 and Xe-NH3 we 
develop full potential energy surfaces, while for the other systems we characterize the stationary 
points (global and local minima and saddle points) . We also calculate isotropic and anisotropic Van 
der Waals Ce coefficients for all the systems. The potential energy surfaces for interaction of NH3 
with alkali-metal and alkaline-earth atoms all show deep potential wells and strong anisotropies. 
The well depths vary from 887 cm" 1 for Mg-NH 3 to 5104 cm" 1 for LI-NH3. This suggests that all 
these systems will exhibit strong inelasticity whenever inelastic collisions are energetically allowed 
and that sympathetic cooling will work only when both the atoms and the molecules are already in 
their lowest internal states. Xe-NH3 is more weakly bound and less anisotropic. 



I. INTRODUCTION 

There is great interest at present in producing samples 
of cold molecules (below 1 K) and ultracold molecules 
(below 1 mK). Such molecules have many potential ap- 
plications. High-precision measurements on ultracold 
molecules might be used to measure quantities of fun- 
damental physics interest, such as the electric dipolc mo- 
ment of the electron [l[ and the time-dependence of fun- 
damental constants such as the electron/proton mass ra- 
tio 0. Ultracold molecules are a stepping stone to ul- 
tracold quantum gases Q and might have applications in 
quantum information and quantum computing Q. 

There are two basic approaches to producing ultra- 
cold molecules. In direct methods such as Stark decel- 
eration 0, 0] and helium buffer-gas cooling 0] , preexist- 
ing molecules are cooled from higher temperatures and 
trapped in electrostatic or magnetic traps. In indirect 
methods Q , laser-cooled atoms that are already ultracold 
are paired up to form molecules by either photoassocia- 
tion Q or tuning through magnetic Feshbach resonances 

Indirect methods have already been used extensively 
to produce ultracold molecules at temperatures below 1 
/jK. However, they are limited to molecules formed from 
atoms that can themselves be cooled to such tempera- 
tures. Direct methods are far more general than indirect 
methods, and can in principle be applied to a very wide 
range of molecules. However, at present direct meth- 
ods are limited to temperatures in the range 10-100 mK, 
which is outside the ultracold regime. There is much 
current research directed at finding second-stage cooling 
methods to bridge the gap and eventually allow directly 
cooled molecules to reach the region below 1 jjK where 
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quantum gases can form. 

One of the most promising second-stage cooling meth- 
ods that has been proposed is sympathetic cooling. The 
hope is that, if a sample of cold molecules in brought into 
contact with a gas of ultracold atoms, thermalization will 
occur and the molecules will be cooled towards the tem- 
perature of the atoms. Sympathetic cooling has already 
been used successfully to cool atomic species such as 6 Li 
PH and 41 K [l2j but has not yet been applied to neutral 
molecules. 

Sympathetic cooling relies on thermalization occurring 
before molecules are lost from the trap. Thermalization 
requires elastic collisions between atoms and molecules to 
redistribute translational energy. However, electrostatic 
and magnetic traps rely on Stark and Zccman splittings 
and trapped atoms and molecules are not usually in their 
absolute ground state in the applied field. Any inelastic 
collision that converts internal energy into translational 
energy is likely to kick both colliding species out of the 
trap. The ratio of elastic to inelastic cross sections is 
thus crucial, and a commonly stated rule of thumb is 
that sympathetic cooling will not work unless elastic cross 
sections are a factor of 10 to 100 greater than inelastic 
cross sections for the states concerned. 

Inelastic cross sections for atom-atom collisions are 
sometimes strongly suppressed by angular momentum 
constraints. In particular, for s-wave collisions (end-over- 
end angular momentum L = 0), pairs of atoms in spin- 
stretched states (with the maximum possible values of the 
total angular momentum F and its projection |Mf|) can 
undergo inelastic collisions only by changing L. Cross 
sections for such processes are very small because, for 
atoms in S states, the only interaction that can change L 
is the weak dipolar coupling between the electron spins. 
However, for molecular collisions the situation is differ- 
ent: the anisotropy of the intcrmolccular potential can 
change L, and this is usually much stronger than spin- 
spin coupling. 

It is thus crucial to investigate the anisotropy of the 
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interaction potential for systems that are candidates 
for sympathetic cooling experiments. In experimen- 
tal terms, the easiest systems to work with are those 
in which molecules that can be cooled by Stark de- 
celeration (such as NH3, OH and NH) interact with 
atoms that can be laser-cooled (such as alkali-metal and 
alkaline-earth atoms). There has been extensive work 
on low-energy collisions of molecules with helium atoms 
[H [3 GJ, El Q3, 03, 0, but relatively little on colli- 
sions with alkali-metal and alkaline-earth atoms. Soldan 
and Hutson [2(| investigated the potential energy sur- 
faces for Rb + NH and identified deeply bound ion-pair 
states as well as weakly bound covalent states. They 
suggested that the ion-pair states might hinder sympa- 
thetic cooling. Lara et al. [2ll . I22I ] subsequently calcu- 
lated full potential energy surfaces for Rb + OH. for both 
ion-pair states and covalent states, and used them to in- 
vestigate low-energy elastic and inelastic cross sections, 
including spin-orbit coupling and nuclear spin splittings. 
They found that even for the covalent states the potential 
energy surfaces had anisotropics of the order of 500 cm -1 
and that this was sufficient to make the inelastic cross sec- 
tions larger than inelastic cross sections at temperatures 
below 10 mK. Tacconi et al. [23| have recently carried out 
analogous calculations on Rb + NH, though without con- 
sidering nuclear spin. There has also been a considerable 
amount of work on collisions between alkali metal atoms 
and the corresponding dimers 0, [H, ES HE H, HI . 

One way around the problem of inelastic collisions is 
to work with atoms and molecules that are in their abso- 
lute ground state in the trapping field. However, this is 
quite limiting: only optical dipole traps and alternating 
current traps [30j ] can trap such molecules. It is therefore 
highly desirable to seek systems in which the potential 
energy surface is only weakly anisotropic. The purpose 
of the present paper is to survey the possibilities for col- 
lision partners to use in sympathetic cooling of NH3 (or 
ND3), which is one of the easiest molecules for Stark de- 
celeration. 

Even if sympathetic cooling proves to be impractical 
for a particular system, the combination of laser cool- 
ing for atoms and Stark deceleration for molecules offers 
opportunities for studying molecular collisions in a new 
low-energy regime. For example, experiments are under 
way at the University of Colorado [31| to study collisions 
between decelerated NH3 molecules and laser-cooled Rb 
atoms. 

There alkali-metal atom + NH3 systems have not been 
extensively studied theoretically, though there has been 
experimental interest in the spectroscopy of Li-NH3 com- 
plex as a prototype metal atom-Lewis base complex [32J . 
Lim et al. (33j ] recently calculated electrical properties 
and infrared spectra for complexes of NH3 with alkali- 
metal atoms from K to Fr and gave the equilibrium struc- 
tures of their global minima. However, to our knowledge, 
no complete potential energy surfaces have been pub- 
lished for any of these systems. The alkaline-earth + 
NH3 have been studied even less, and except for an early 



study of the Be-NH3 system [34[ there are no previous 
results available. 



II. 



AB INITIO METHODS 



The interaction energy of two monomers A and B is 
defined as 
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where E^f is the total energy of the dimcr and E^ ot 
and Efo t are the total energies of the isolated monomers. 
Since the interaction energy is dominated at long range 
by intermolecular correlation (dispersion), ah initio cal- 
culations of the interaction energy must include elec- 
tronic correlation effects at the highest possible level [35[ 
and must be carried out with large basis sets augmented 
by diffuse functions. At present, the coupled-cluster (CC) 
method with single, double and nonitcrativc triple excita- 
tions (CCSD(T)) provides the best compromise between 
high accuracy and computational cost. In the present pa- 
per, we carry out coupled-cluster calculations using the 
Molpro package J3f|. All interaction energies are cor- 
rected for basis-set superposition error (BS SE) with the 
counterpoise method of Boys and Bcrnardi [3 71 ] . 

Standard coupled-cluster methods are reliable only 
when the wavefunction is dominated by a single elec- 
tronic configuration This is often an issue for molecu- 
lar systems with low-lying excited states. In order to 
check the reliability of CC calculations, it is necessary to 
monitor the norm of T\ operator (38j (measured by the 
Tl diagnostic). In the case of metal-NH3 systems this 
is relatively large, especially when the atom approaches 
the lone pair of the NH3 molecule, but the convergence of 
the CC equations is fast and the converged CCSD results 
are very close to benchmark multircfcrcnce configuration 
interaction (MRCI-SD) calculations with size-extensivity 
corrections. Thus we consider the CC results reliable. 

To understand the origin of the intermolecular forces 
we also consider the interaction energies obtained at the 
Hartree-Fock level, which neglects electron correlation 
and thus provides information about the role of disper- 
sion and other correlation effects. For some systems we 
also analyze the components of the intermolecular in- 
teractions using symmetry-adapted perturbation theory 
H| (SAPT). The first-order SAPT corrections (electro- 
static and exchange terms) are computed at the Hartree- 
Fock level, while the dispersion energy is evaluated in the 
coupled Hartree-Fock approximation . These calcula- 



tions are carried out using the SAPT2006[41| program. 

We are interested principally in the collisions of cold 
ammonia molecules with atoms at energies that are much 
too low for vibrational excitation to occur. Such colli- 
sions are governed by an effective potential that is vi- 
brationally averaged over the ground-state vibrational 
wavefunction of NH3 . For the present purpose it is ade- 
quate to represent this by a potential calculated with the 
NH3 molecule frozen at a geometry that represents the 
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ground state. In the present paper we use a geometry de- 
rived from the high-resolution infrared spectra [4^ |: the 
molecule is taken to have C% v symmetry with N-H bond 
lengths of 1.913 a and an H-N-H angle of 106.7°. Inter- 
molecular geometries are specified in Jacobi coordinates: 
R is the distance from the center of mass of NH3 to the 
atom, while 8 is the angle between the intermolecular 
vector and the C3 axis of the NH3 molecule (with 8 = 0° 
corresponding to the atom approaching towards the lone 
pair of NH3). Finally, \ is the dihedral angle between the 
plane containing the C3 axis and an NH bond and that 
containing the C3 axis and the intermolecular vector. 

Table U gives the lowest excitation energies, dipole po- 
larizabilities and ionization energies of the atoms studied 
in this paper. The neutral alkali-metal and alkaline-earth 
atoms (denoted below as A and Ae, respectively) have 
particularly low excitation energies, resulting from small 
separations between energy levels corresponding to ns 
and np or (n — l)d configurations. Since the gap between 
the ground and excited states is small, the atoms have 
very large polarizabilities. Hence, we expect particularly 
strong induction and dispersion interactions. The alkali- 
metal and alkaline-earth atoms also have low ionization 
energies E\ . Since the atomic orbital wavefunctions van- 

1 11 

ish at long range as exp(— E i r), the wavefunctions and 
densities are very diffuse, and this causes large overlap 
between monomers even at relatively large separations. 
Finally, because of the low ionization energies, alkali- 
metal and alkaline-earth atoms have a strong tendency 
to form charge-transfer complexes. 

The basis sets used in the ab initio calculations are 
as follows. For Be, Li, Mg, Na, Ca atoms we use all- 
electron cc-pVTZ basis sets augmented by even-tempered 
diffuse exponents, while for potassium we use the CVTZ 
basis set of Feller JH^]. For Rb, Sr and Xe we han- 
dle only the outermost electrons explicitly, with the core 
electrons represented by effective core potentials (ECPs). 
For Rb we use the small-core effective core potential 
ECP28MWB with a basis set based on that of Ref. [H, 
which was optimized to recover the static dipole polar- 
izability. We modified this slightly to account better for 
intramonomer electronic correlation effects by removing 
0.07 / and adding 0.001049 s, 0.0024 p, 4.5,0.016667 d, 
1.9,0.655 / and 0.95,0.3167 g functions. The basis set 
for Sr is taken from Ref. l54l . For Xe we use the basis set 
given by Lozeille et al. [50j, which was found to be ex- 
cellent for polarizabilities and hyperpolarizabilities. For 
each system we added a set of midbond functions with 
exponents sp: 0.9,0.3,0.1, df 0.6,0.2 to improve the rep- 
resentation of the dispersion energy in the region of the 
Van der Waals minimum. 



III. RESULTS AND DISCUSSION 

The potential energy surface for an atom-NBi3 system 
is a function of the intermolecular distance R and two 
angles 8 and X - However, functions of 3 variables are 



difficult to represent graphically. It is convenient to rep- 
resent the x-dependence in the form 

00 

V(R,6, X ) = Y, V M8)™s3k X . (2) 

fc=0 

To reduce the computational effort we calculate the in- 
teraction potential only for \ = 0° and x = 60° and 
approximate the leading terms Vo(R,8) and Vs(R,8) by 
sum and difference potentials, 

V (R,9) = ±[V(R,9,0)+V(R,8,60°)] 
V 3 (R,6) = ±[V(R,6,0)-V(R,8,60°)\. (3) 

Vo can be viewed as the interaction potential averaged 
over Xi while V3 describes the leading anisotropy of the 
potential with respect to rotation about the C3 axis of 
NH 3 . 

A. Alkali-metal atom + NH3 interactions 

The potential energy surface for Rb-NH3 is shown in 
Figure [TJ CCSD(T) calculations were carried out at 
\ = 0° and 60° , at values of 8 corresponding to a 20-point 
Gauss-Lobatto quadrature. The grid included R values 
from 3.5 to 12 ao with a step of 0.5 do, and from 12 to 
15 ao with the step of 1 ao. There is a deep minimum 
(1862 cm -1 ) at R = 5.90 ao and 8 = 0, corresponding to 
approach of Rb towards the NH3 lone pair. The poten- 
tial is much shallower at other geometries, with a saddle 
point near 8 = 110° and a shallow secondary minimum 
at 8 = 180°. The anisotropy with respect to rotation of 
NH3 about the C3 axis (x) is relatively weak, at least 
in the low-energy classically allowed region defined by 
V (R,6)<0. 

The overall shape of the other A-NH3 potentials is 
quite similar. In each case there is a deep minimum 
around 8 = 0° and a shallow secondary minimum for 
8 = 180°. Table Ull gives the well depths and equilibrium 
distances. For the alkali metals the well depth of the 
global minimum decreases down the periodic table, from 
5104 cm -1 for Li to 1862 cm -1 for Rb, and the equilib- 
rium distance increases from 3.91 ao for Li to 5.90 ao for 
Rb. The changes in the properties of the shallow sec- 
ondary minima are much smaller, with well depths close 
to 100 cm -1 for all the alkali metals. Our results for 
the species containing K and Rb are in good agreement 
with the CCSD(T) calculations of Lim et al. [33]; they 
obtained slightly different values of the binding energies 
of K-NH 3 and Rb-NH 3 (2210 cm" 1 and 1950 cm" 1 , re- 
spectively), but their results are not corrected for BSSE. 
It should also be noted that their binding energies are for 
relaxed NH 3 geometries. 

The deep wells and large anisotropies of the A-NH3 
potentials will produce strong coupling between the dif- 
ferent NH3 rotational states during collisions. All these 
systems are therefore likely to have large inelastic cross 
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TABLE I: Properties of alkali-metal, alkaline-earth and Xe atoms important to interaction potentials. Note that for alkali-metal 
atoms the lowest excitation energy corresponds to 2 S — > 2 P\/2 excitation and for alkaline-earth and Xe atoms to 1 S — > 3 Po- 
The excitation and ionization energies are take from NIST Handbook of Basic Atomic Spectroscopic Data [43j. 



Atom 


dipole polarizability 
(ag) 


Ref. 


Ce coefficient 
(E h a 6 ) 


Ref. 


lowest excitation energy 
(cm" 1 ) 


ionization energy 
(cm" 1 ) 


Li 


164 


44 


1395 


45 


14904 


43487 


Na 


162 


46 


1561 


45 


16956 


41449 


K 


293 


46 


3921 


47 


12985 


35010 


Rb 


319 


46 


4707 


48 


12578 


33691 


Be 


37.7 


45 


213 


45 


21978 


75193 


Mg 


71 


45 


629 


45 


21850 


61671 


Ca 


159 


45 


2221 


49 


15157 


49305 


Sr 


200 


45 


3250 


45 


14317 


45932 


Xe 


27.3 


50 


286 


51 


67068 


97834 
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FIG. 1: The interaction potential of Rb-NH 3 from CCSD(T) 
calculations: Vo(R,6) component (upper panel) and Vs(R, 6) 
component (lower panel). Contours are labelled in cm -1 . To 
aid visualization, V3 is plotted only in the energetically acces- 
sible region defined by Vo < 0. 



sections. It is thus unlikely that sympathetic cooling 
of NH3 with alkali-metal atoms will be successful unless 
both the atoms and the molecules are already are in their 
lowest internal states. 



TABLE II: Equilibrium distances and well depths for alkali- 
metal atom + NH3 systems from CCSD(T) calculations. 



6 = 0° 9 = 180 





Re (do) 


D e (cm" 1 ) 


-Re (ao) 


D e (cm" 1 ) 


Li 


3.91 


5104 


7.86 


104.8 


Na 


4.73 


2359 


8.33 


98.2 


K 


5.52 


2161 


8.90 


99.6 


Rb 


5.90 


1862 


8.89 


110.2 



B. Alkaline-earth atom + NH3 interactions 



We originally hoped that the potentials for systems 
containing alkaline-earth atoms would be more weakly 
bound and less anisotropic than for those containing 
alkali-metal atoms . However, this proved not to be the 
case, at least for the heavier alkaline-earth atoms that are 
most suitable for laser cooling. The results for the Ae- 
NH3 systems are summarized in Table IIIII The shapes 
of the potential energy surfaces are generally similar to 
those for A-NH3 systems. For Ca and Sr, the depths 
of the global minima are 3229 and 3141 cm^ 1 respec- 
tively; these are both deeper than for the corresponding 
alkali-metal atom. For Mg, however, the well depth is 
considerably shallower at only 887.5 cm -1 . The min- 
ima corresponding to approach at the hydrogen end of 
NH3 arc slightly deeper than for the alkali metals, rang- 
ing from 115.7 for Mg to 131.6 cm -1 for Sr. On the 
other hand, the interaction potential for Be-NH3 resem- 
bles those for Ca-NH3 and Sr-NH3 more than that for 
Mg-NH3: the global minimum is 1973 cm -1 deep, while 
the dispersion-bound minimum is 100.5 cm -1 deep. The 
equilibrium distance for Be-NH3 at 8 = 0° (3.57 ao) is 
also much shorter than for the other AC-NH3 systems, 
and is comparable to that for Li-NH3. 
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TABLE III: Equilibrium distances and well depths for 
alkaline-earth atom + NH3 systems from CCSD(T) calcula- 
tions. 




e = o° 


e 


= 180° 




R e (a ) D e (cm' 1 ) 


Re (do) 


D e (cm" 1 ) 


Be 


3.57 1973 


7.61 


100.5 


Mg 


4.83 887.5 


8.20 


115.7 


Ca 


4.92 3229 


8.85 


129.1 


Sr 


5.22 3141 


9.06 


131.6 


TABLE IV: The interaction energies (in cm -1 ) for L1-NH3, 
Ca-NH3 and Mg-NH3 at different levels of electronic correla- 
tion, for geometries corresponding to the global and secondary 
minima. 






CM 






HF 


CCSD 


CCSD(T) 


Li 


-4405 


-5022 


-5104 


Ca 


-2152 


-2937 


-3229 


Mg 


260 


-590 


-888 






LM 






HF 


CCSD 


CCSD(T) 


Li 


248 


-54 


-105 


Ca 


244 


-47 


-129 


Mg 


155 


-54 


-116 



C. Origin of bonding in metal-atom + NH3 

systems 

It is important to understand the large difference be- 
tween the metal-lone pair bond energies between Mg and 
the other Group 1 and 2 atoms considered here. Table ITVl 
gives the interaction energies in the global and secondary 
minima at the Hartree-Fock, CCSD and CCSD(T) levels 
for Li-NH 3 , Mg-NH 3 and Ca-NH 3 . For all these systems 
the Hartree-Fock interaction energies are positive for the 
shallow secondary minima, indicating that the shallow 
wells are dominated by dispersion forces. At the global 
minima, however, Mg-NH3 is repulsive at the Hartree- 
Fock level while the other two systems are strongly at- 
tractive. There is thus strong chemical bonding in Li- 
NH 3 and Ca-NH 3 that is absent in Mg-NH 3 . 

The qualitative differences between Mg and the other 
atoms can be understood if we consider how the energy 
of the highest occupied molecular orbital (HOMO) differs 
for the different atom-NH 3 systems. Fig.[2]shows the two 
highest occupied molecular orbitals of each system. As 
we separate the monomers to infinity, these two orbitals 
became HOMOs of the atom and the NH 3 molecule. For 
any alkali-metal atom, the strong A-NH 3 bond can be 
explained in terms of LCAO-MO theory as a chemical 
bond of order one half, since we have a doubly occupied 
bonding orbital and a singly occupied antibonding or- 



FIG. 2: The pattern of molecular orbitals for a) Li-NH3, b) 
Ca-NH 3 , c) Mg-NH 3 near their global minima. The HOMOs 
of NH3 and of the metal atoms form bonding and antibonding 
orbitals. Note the small change in the HOMO energy for the 
LLNH3 and Ca-NH3 systems and the much larger change for 
Mg-NH 3 . 



bital [see Fig. [2] a)]. However, this explanation does not 
apply to the alkaline-earth atoms, where the antibonding 
orbital is doubly occupied. The net bonding in Ca-NH 3 
arises because the bonding orbital is shifted down in en- 
ergy considerably more than the antibonding orbital is 
shifted up. Conversely, in Mg-NH 3 , the contributions 
from the bonding and antibonding orbitals are closely 
balanced. The difference can probably be attributed to 
the participation of np orbitals; as shown in Table U the 
S — > P splitting is considerably smaller in Ca than in 
Mg. Thus Mg-NH 3 is bound mainly by dispersion forces 
whereas Ca-NH 3 has substantial chemical bonding. 

Different considerations apply to the Be atom, which is 
a notoriously difficult case for electronic structure theory 
[55^ Although the potential energy surfaces are qual- 
itatively similar for the Be-NH 3 , Ca-NH 3 and Sr-NH 3 
systems at the CCSD(T) level, the origin of the strong 
bonding is probably different in Be-NH 3 . In this case 
the Hartree-Fock and CCSD potential energy curves for 
9 = 0° show a double-minimum structure, with a shal- 
low long-range minimum separated from the global min- 
imum by a barrier. This suggests a sudden change in 
chemical character as the Be atom approaches N. At the 
Hartree-Fock level the maximum has an energy of 730 
cm -1 at R = 4.88 clq. The long-range minimum at the 
Hartree-Fock level is 18.4 cm -1 deep at R = 9.02 do, 
while at the CCSD level it is 138 cm -1 deep at R = 6.5 
ao. Despite this peculiar behavior, the CC calculations 
showed no convergence problems or unusually large Tl 
diagnostics. However, our results for Be-NH 3 disagree 
with those of Chalasiriski and coworkers [34[ , who carried 
out fourth-order Moller-Plesset (MP4) calculations and 
found a global minimum that corresponds to the outer 
minimum on the CCSD potential energy curve. They did 
not find the inner minimum, which turned out to be the 
global minimum in our calculations. 

As mentioned before, the feature of the potential en- 
ergy surfaces that is important for elastic/inelastic col- 
lision ratios is the anisotropy. In order to understand 
the origin of the anisotropics better, we carried out ad- 
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FIG. 3: Electrostatic energy of Na (squares) and Mg (circles) 
atoms interacting with NH3 as a function of 9 for 7? = 600. 
The energy is averaged over %■ 



FIG. 5: Dispersion energy of Na (squares) and Mg (circles) 
atoms interacting with NH3 as a function of 9 for R = 600. 
The energy is averaged over \- 
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FIG. 4: First-order exchange energy of Na (squares) and Mg 
(circles) atoms interacting with NH3 as a function of 9 for 
R — 6ao. The energy is averaged over \. 



ditional calculations based on symmetry-adapted pertur- 
bation theory (SAPT). Figures [3J 2J and[5]sliow the elec- 
trostatic, first-order exchange and dispersion components 
of the interaction energy Vq for Na-NH3 and Mg-NH3, 
averaged over x as m Eq. (J3j) . The calculations were 
performed at a fixed R value of 6 ao, which is in an 
attractive region for 9 = 0° and a repulsive region for 
9 = 180°. Figs. EJ SI and [5] show clearly that it is the 
first-order interaction energy that is responsible for most 
of the anisotropy in the valence overlap region. This is 
caused by a very large difference between the electrostatic 
attraction near the lone-pair site and near hydrogen sites 
(see Fig. [3J . This difference is significantly larger than 
that in the exchange energy The anisotropy of the dis- 
persion interaction (plotted in Fig. [5]), is even weaker. 

The anisotropy V3 of all three components of the inter- 
action energy with respect to x IS shown for Na-NH3 and 
Mg-NH 3 in Fig. [6] The exchang e energy is very strongly 
anisotropic, especially for Mg-NH3. The large difference 
in exchange energy between Mg-Nfi3 and Na-NH3 can be 
explained by the closed-shell character of the Mg atom 
and the much stronger Pauli repulsion between hydro- 
gens of NH 3 and Mg. The electrostatic and dispersion 
contributions to V3 are much more similar for the two 
systems. 
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FIG. 6: Anisotropy of the electrostatic, exchange and disper- 
sion contributions to the interaction energy, with respect to 
rotation about the C3 axis of NH3, as a function of 9, for 
R = 6 ao, for Na-NH3 (upper panel) and Mg-NH3 (lower 
panel) . 



D. Xe + NH3 interaction 

All the metal-NH3 systems investigated above have dis- 
appointingly large anisotropics. It is likely that all of 
them will exhibit large inelastic cross sections for any 
initial state where inelasticity is possible. We there- 
fore decided to consider other possible collision partners 
for sympathetic cooling of NH3. Barker [|36| has sug- 
gested an experiment in which Xe is first laser-cooled 
in its metastable 3 P2 state and then transferred to its 
ground 1 5'o state by laser excitation followed by sponta- 
neous emission. Since ground-state Xe has a fairly large 
dipole polarizability, it can be held in an optical dipolc 
trap and might be used for sympathetic cooling. In this 
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subsection wc investigate the Xe-NH3 interaction in order 
to evaluate its potential in this respect. 

Interactions between noble gases and ammonia have 
been studied extensively. The interaction between He 
and NH3 is important in understanding th e sp ectroscopy 
of NH3 molecules in helium nanodroplets [57[ . The most 
recent ab initio calculations of Hodges and Wheatley 
[HI, [59| gave a global minimum about 33 cm -1 deep at 
R = 6 a , 9 = 90° and x = 60°. The interaction of Ar 
with NH3 has been studied even more extensively, both 
experimentally fsoj ] and by ab initio methods f6ll [6 p . In- 
version of vibration-rotation-tunnelling spectra [60j gave 
a minimum 147 cm -1 deep at R = 6.5 ao, 9 = 97° 
and x = 60°, while the ab initio MP 4 (fourth-order 
M0ller-Plesset) calculations of Tao and Klemperer [62j 
gave a global minimum 130 cm -1 deep at R = 6.85 ao, 
9 = 90° and \ = 60°. The Ne-NH 3 system was inves- 
tigated through MP4 calculations by van Wijngaarden 
and Jager [63|, who obtained a global minimum 63 cm -1 
deep at R = 6.1 a , 9 = 90° and x = 60°. For Kr- 
NH3, Chalasihski et al. [64| obtained a global minimum 
108 cm- 1 deep at R = 7.2 a , 9 = 100° and x = 60°. 
However their results were based on calculations at the 
MP2 level and may not reproduce the dispersion energy 
accurately. 

Fig.[7]shows the interaction potential for Xe-NHs from 
our CCSD(T) calculations. The potential energy sur- 
face differs qualitatively from those for metal-NH3 po- 
tentials studied in the previous subsection, and behaves 
analogously to those for other Rg-NH3 systems. The Vb 
surface for XC-NH3 has only one minimum, 173.5 cm -1 
deep, at R = 7.65 ao and 9 = 66°. The global mini- 
mum for the non-expanded surface is 196.8 cm _1 deep, 
at R = 7.35 a , 6 = 81° and X = 60°. There are saddle 
points at both C% v geometries. For 9 = the saddle point 
is 166.2 cm- 1 deep at R = 7.73 a , while for 9 = 180° 
the saddle point is 134.1 cm -1 deep at R = 7.93 ao. The 
major binding arises from the dispersion energy, and at 
the Hartree-Fock level we observe only a small attraction 
(a few cm -1 ) at large distances, due to weak induction 
forces which behave asymptotically as —C§R~ 6 . Near 
the Van der Waals minimum predicted by CCSD(T), the 
SCF energy is repulsive. 

The Vb surface for Xe-NH3 system thus has an 
anisotropy of only about 60 cm -1 between the poten- 
tial minimum and the higher of the two saddle points. 
This is considerably smaller than for Rb-OH or any of 
the metal-NH3 systems studied here, but still substan- 
tial compared to the rotational constant of NH3, b = 6.35 
cm -1 for rotation about an axis perpendicular to C3. 



E. Long-range forces 

Long-range forces are very important in cold and ultra- 
cold collisions. We therefore carried out separate calcu- 
lations of the Van der Waals coefficients for the interac- 
tions. The isotropic Cg.o and anisotropic Cg,2 dispersion 




FIG. 7: The interaction potential of Xe-NH 3 from CCSD(T) 
calculations: Vo{R, 9) component (upper panel) and Vs(R, 6) 
component (lower panel). Contours are labelled in cm -1 . To 
aid visualization, V3 is plotted only in the energetically acces- 
sible region defined by Vb < 0. 



coefficients for the interaction of atom A and symmetric 
top molecule B may be written in terms of the dynamic 
polarizabilitics of the monomers, evaluated at imaginary 
frequencies, 



°6,0 



a a (iu)a b (iu)du; 



C, 



disp 



+ 3C- 



(4) 



where a = ^(2a xx + a zz ) is the isotropic polarizability 
and Aa = a zz ~a xx is the polarizability anisotropy. The 
induction contributions to the Van der Waals coefficients 



Cind 
k n 



Cind 



a A fl 



(5) 



where the dipole moment fi is 0.579 eao for NH3 [651 ] . 

The integrals in Eqs. [4] were evaluated using the 
method given by Amos et al. (66| . The dynamic po- 
larizabilitics of NH3 were obtained using coupled Kohn- 
Sham theory with the asymptotically corrected PBE0 
functional [67| and d-aug-cc-pVTZ basis sets. To get 
the dynamic polarizabilitics for the alkali- metal atoms, 
we adjusted the fraction of exchange, exact exchange and 
correlation fraction in the PBE0 functional in such a way 
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TABLE V: Van der Waals dispersion and induction coeffi- 
cients for A-NH3 and Ae-NH3 systems. All values are in 
atomic units, E^a^. 









QCo — Cm 


Li 


224 


7.2 


55.0 


Na 


258 


7.4 


54.3 


K 


378 


11.6 


98.2 


Rb 


416 


12.5 


106.9 


Be 


121 


2.3 


12.7 


Mg 


200 


4.4 


24.0 


Ca 


342 


8.4 


53.4 


Sr 


413 


10.2 


67.4 


Xe 


161 


0.94 


9.1 



as to recover the atom-atom Cq coefficients (see Table HJ. 
Our coupled Kohn-Sham program does not allow us to 
use core potentials to calculate dynamic polarizabilitics. 
For Rb we therefore performed all-electron calculations 
with the pVTZ basis set of Sadle.j [68j combined with the 
Douglas-Kroll approximation [69| . The dynamic polariz- 
abilitics obtained in this way were tested by comparing 
Cq coefficients for A-Ar and A-Xe systems with those ob- 
tained by Mitroy and Zhang [7(i| . The maximum error 
was found to be +6.3% (for Na-Xe) while the average 
error is less than +3%. For alkaline-earth and Xc atoms 
the frequency-dependent dipole polarizabilities were ob- 
tained from time-independent coupled-cluster linear re- 
sponse functions [7l|, [72| . 

The resulting Cq coefficients are shown in Table fVl For 
the alkali- metal and alkaline-earth atoms, the isotropic 
dispersion coefficients Cg q P are fairly large because of 
the large atomic polarizabilitics. The anisotropics in the 
dispersion coefficients are much smaller, because of the 
small polarizability anisotropy of NH3 (2.1 <2q) compared 
to its isotropic polarizability (14.6 Oq). The induction 
Van der Waals coefficients are large, and account for 10- 
25% of the total C 6j0 and 70-90% of the total C 6 , 2 . It 
may be noted that C 6 q P for Rb-NH 3 is somewhat larger 

than CgQ P for Rb-OH [22|. As one might expect, the Xe- 
NH3 long-range interaction has slightly different charac- 



ter from the A- and Ae-NH3 systems. The C 6 q P coeffi- 
cient is still large, but the total anisotropy (in particular 
the dispersion anisotropy) is much smaller. 



IV. CONCLUSIONS 

We have investigated the intcrmolccular potential en- 
ergy surfaces for interaction of NH3 with several different 
atoms that might be used for sympathetic cooling. For 
interaction with all the alkali-metal and alkaline-earth 
atoms, we found deep minima and strong anisotropies. 
The shallowest potential is for Mg-NH3, but even there 
the anisotropy in the well depth is close to 800 cm -1 . 
This is likely to cause strong inelastic collisions for all 
initial states for which they are energetically allowed. Ac- 
cordingly, wc consider that none of the alkali metals and 
alkaline earths arc good prospects for sympathetic cool- 
ing of NH3 unless both the atoms and the molecules are 
in their lowest states in the trapping field. This suggests 
that sympathetic cooling would need to be carried out in 
either optical or alternating current traps. 

A somewhat more promising system for sympathetic 
cooling is Xe-NH3, for which the global minimum is cal- 
culated to be 196.8 cm" 1 deep at an off-axis geometry. 
The XC-NH3 system is relatively weakly anisotropic, with 
the saddle points for C% v geometries only 30.6 and 62.7 
cm" 1 higher than the global minimum. In future work we 
will use the interaction potential to calculate low-energy 
elastic and inelastic cross sections, in order to predict 
whether sympathetic cooling of NH3 by Xe is likely to be 
feasible. 

Even if sympathetic cooling proves to be impossible for 
these systems, there is much to be learnt from collisions 
between velocity-controlled beams of molecules and laser- 
cooled atoms. There are opportunities to explore low- 
energy inelastic processes in novel collisional regimes and 
to probe scattering resonances in unprecedented detail. 
We therefore intend to use the potential energy surfaces 
developed here to carry out inelastic collision calculations 
to explore these effects and assist in the interpretation of 
collision experiments. 
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